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We find novel phase transitions and critical phenomena that occur only outside the linear-response regime of 
current-driven nonequilibrium states. We consider the strongly-interacting (3 + 1) -dimensional M = 4 large- 
N c SU(Nc) supersymmetric Yang-Mills theory with a single flavor of fundamental TV — 2 hypermultiplet as 
a microscopic theory. We compute its nonlinear non-ballistic quark-charge conductivity by using the AdS/CFT 
correspondence. We find that the system exhibits a novel nonequilibrium first-order phase transition where the 
conductivity jumps and the sign of the differential conductivity flips at finite current density. A nonequilibrium 
critical point is discovered at the end point of the first-order regime. We propose a nonequilibrium steady-state 
analogue of thermodynamic potential in terms of the gravity-dual theory in order to define the transition point. 
Nonequilibrium analogues of critical exponents are proposed as well. The critical behavior of the conductiv- 
ity is numerically confirmed on the basis of these proposals. The present work provides a new example of 
nonequilibrium phase transitions and nonequilibrium critical points. 
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Introduction. — Nonequilibrium physics is one of the cen- 
tral subjects in modern physics. Although the linear response 
theory provides a computational framework of the transport 
coefficients at the vicinity of thermal equilibrium, its exten- 
sion to the nonlinear regime is still a great challenge. A 
key question in nonequilibrium physics is how to extract the 
macroscopic physics from the underlying microscopic theo- 
ries when the systems are out of equilibrium. 

Recently, a great deal of attention has been paid to the 
anti-de Sitter space (AdS)/conformal field theory (CFT) corre- 
spondence fUll]. AdS/CFT correspondence is a map between 
a strongly interacting quantum gauge theory and a classical 
gravity. The map is essentialy established at the level of mi- 
croscopic theory. However, it also provides a new picture for 
statistical mechanics. The macroscopic physics, such as ther- 
modynamics, of the gauge particles appears naturally in terms 
of the black hole physics on the gravidy side The point is 
that the coarse graining of the gauge theory is accomplished 
in the gravity side just by solving the classical equations of 
motion. This simplification tempts us to apply AdS/CFT to 
many-body physics out of equilibrium. 

Usually, a drawback of AdS/CFT is that we need to deal 
with an idealized gauge theory that is not exactly realized in 
nature at the microscopic level, in order to ensure that the cor- 
respondence is well defined. However, the aim of nonequilib- 
rium statistical mechanics is to describe macroscopic proper- 
ties common to a wide range of many -body systems regardless 
of the details of each microscopic theory. Our ultimate goal is 
to obtain new information on such macroscopic physics that 
is shared by the actual systems in nature. Phase transition and 
critical phenomena are ideal places to look for such informa- 
tion in light of their possible universality. 

In this Letter, we study nonequilibrium phase transitions |01 
by using the AdS/CFT correspondence. Our system consists 
of a strongly interacting gauge-theory plasma driven to the 
nonequilibrium steady state (NESS) by a constant current. 
We discover novel nonequilibrium phase transitions and a 



nonequilibrium critical point in the nonlinear regime. We de- 
velop a formalism to analyze the nonequilibrium phase transi- 
tions as well: we propose NESS analogues of thermodynamic 
potential and critical exponents. 

Our phase transitions are associated with the nonlinear 
conductivity of our system. It has been shown in Ref. 
that the gauge-theory plasma, which will be defined in detail 
later, exhibits negative differential conductivity (NDC) in low- 
current-density regions. This NDC is categorized as S-shaped 
NDC (SNDC) in Ref. GJ. Here, the differential conductivity 
is defined as dJ/dE, where J and E are the current density 
and the external electric field acting on the charge carriers, 
respectively. The difference from the conventional conduc- 
tivity a = J/E is that the differential conductivity here can 
be either negative or positive, whereas a cannot be negative. 
In fact, SNDC has been experimentally observed in various 
current-driven systems Q, including in systems of strongly 
correlated electrons (for example, see Ref. [0]). It has been 
shown in Ref. 15[] that the NDC is converted to positive differ- 
ential conductivity (PDC) in the high current-density regions. 
The transition from NDC to PDC was observed to be smooth 
in Ref. JSfl: it is a crossover 

In this Letter, we discover first-order and second-order 
phase transitions between NDC and PDC. We also observe 
critical phenomena at the critical point. As far as the au- 
thor knows, the current-driven nonequilibrium phase transi- 
tions between NDC and PDC, and the associated critical point, 
have been reported neither experimentally nor theoretically so 
far: the present work provides a new model of nonequilibrium 
phase transitions and nonequilibrium critical points. 

Microscopic theory. — We choose our microscopic theory 
by asking how much of its AdS/CFT duality is established, 
rather than asking whether the microscopic details are real- 
ized in nature, since we are interested in only the macroscopic 
properties of NESS that may be independent of the micro- 
scopic details. One of the gauge theories whose AdS/CFT 
duality is well-established is the (3 + 1) -dimensional strongly- 
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coupled SU(N C ) Af = 4 super-symmetric Yang-Mills theory 
(TV = 4 SYM) at the large-/Y c limit with a single flavor of 
fundamental Af = 2 hypermultiplet. This is a supersymmet- 
ric cousin of quantum chromodynamics (QCD), but its super- 
symmetry is broken at finite temperatures. The Af = 4 SYM 
sector (gluon sector) contains the gauge particles in the ad- 
joint representation, which we call "gluons" in this Letter. The 
AT = 2 hypermultiplet sector (quark sector) contains particles 
in the fundamental representation and anti-fundamental rep- 
resentation; we call these particles "quarks" and "antiquarks," 
respectively. The quark (antiquark) carries a unit of positive 
(negative) quark charge. These particles play the role of, for 
example, electrons and holes in condensed matter, that is, the 
role of charge carriers. The interaction among these particles 
is mediated by the gluons; the gluons play the role of phonons, 
for example. 

We consider the conductivity associated with the quark 
charge, and it is the quark-charge current that drives the quark 
sector out of equilibrium. NESS is realized in the following 
manner |8[]. We set the large- N c limit that makes the degree 
of freedom (DOE) of the gluon sector [which is 0(N 2 )] suf- 
ficiently larger than that of the quark sector [which is 0(N C )]. 
We also set the gluon sector equilibrium at a definite tempera- 
ture T. The interaction between the gluon sector and the quark 
sector generates a dissipation in the presence of the quark cur- 
rent: the gluons absorb the momentum and the energy of the 
charge carriers. Because of the large DOE, the heat capacity 
of the gluon sector is sufficiently large, and the temperature 
of the gluon sector is well approximated as a constant. The 
gluon sector plays the role of heat bath, and the NESS of the 
quark sector is realized by putting the dissipation into the heat 
bath and the work of the external field in balance. 

Let us specify the conditions we impose on the system. We 
consider a neutral system where the total quark-charge den- 
sity is zero. This means that the finite current is realized by 
equal numbers of quarks and antiquarks flowing in opposite 
directions. The system works as an insulator when E and T 
are sufficiently small compared to the mass of the quark, but 
a strong enough electric field will break the insulation [0 [Toll . 
Our charge carriers are those pair created by the external field 
1 1 111 in the insulation breaking. We assume that the system is 
steady and homogeneous. We also assume that the system has 
an infinitely large volume: note that we are not dealing with 
mesoscopic systems. We consider current-driven phenomena, 
and choose J (the quark-charge current density) as our con- 
trol parameter. In this sense, E, the external field acting on 
the quark charge, is taken as a function of J. We choose the 
electric field (and hence the current) to be in the x direction. 
We employ the natural units c = H = &b = L 

Nonlinear conductivity in AdS/CFT. — The gravity dual of 
our microscopic theory is the so-called D3-D7 system 1112 1. 
The computational technique of nonlinear conductivity has 
been proposed in 11 ill in the framework of the AdS/CFT corre- 
spondence, and we follow it. We sketch the idea of 1 1 lh be low 
to define our notation, and one may consult Refs. Ill llll2ll for 
more details. Our proposals for thermodynamic potential and 



critical exponents shall be given later. 

The gluon sector at finite temperature is mapped to the grav- 
ity theory on a curved geometry, which is a direct product of 
a five-dimensional AdS-Schwarzschild black hole (AdS-BH) 
and S 5 H ■ The metric of the AdS-BH part is given by (B ] 



i (i-*V4) 2 

z 2 l + z 4 /z% 



ds A = - 



■dt z 



i + z 4 /4 



dx 2 



df_ 



(1) 



where z (0 < z < zh) is the radial coordinate of the black 
hole geometry. The AdS-BH has a horizon at z = Zh and a 
boundary at z — 0. The Hawking temperature, which corre- 
sponds to the temperature of the gluon sector (hence that of 
the heat bath), is given by T = \/2/ (ttzh)', x and t denote 
the (3 + l)-dimensional spacetime coordinates on which the 
gauge theory is defined. The metric of the S 5 part is given by 



dnj = de 2 



sin 2 9d^ 2 + cos 2 9dfl 



3! 



(2) 



where < 9 < tt/2, and dfld is the volume element of the 
unit d-dimensional sphere. The radius of the S 5 part is taken 
to be 1 for simplicity. This is equivalent to choosing the 't 
Hooft coupling A of the gauge theory at A = (2_7r) 2 /2. 

The quark sector is mapped to a D7-brane 1411 embedded 
to the above geometry. The three-dimensional part of the D7- 
brane is wrapped on the S 3 part in S 5 , the radius of which is 
given by cos 9: the c onfig uration of the D7-brane is specified 
by the function 6{z) IU2I1 . We employ the probe approxima- 
tion where the backreaction of the D7-brane to the AdS-BH is 
neglected. This is justified at the large-A^ limit, and is con- 
sistent with the picture of the AdS-BH as the heat bath. The 
behavior of 9(z) at the vicinity of the boundary is related to 
the current quark mass m q 1 1 211 



as 



6{z) = m q z + - [(qq)/N + ml/3] z 3 + 0(z 5 ). (3) 

Here, (qq) denotes the chiral condensate |Q . Throughout the 
analysis, we fix m q at a designed value. There is a U(l) gauge 
field on the D7-brane. E and J are related to A x through 

A x {z,t) = -Et + const. + J(2/Y)- 1 z 2 + 0(z 4 ), (4) 

where N = N c /(2tt) 2 in our convention and we have em- 
ployed the gauge d x A t = Ml in . The dynamics of 9(z) 
and A x (z,t) is governed by the D7-brane action, Sd7 = 



J dtd 3 xdz£r>7, where £u7 is explicitly written as JIT 



Ng xx cos 3 9^J \gtt\g xx g zz - 9zz(A x ) 2 + \g tt \{A x ) 2 .(5) 



Here, the prime (the dot) denotes the differentiation with re- 
spect to z (t); gtt, Six and g zz are the components of the in- 
duced metric on the D7-brane, and they are equal to those of 
the metric of AdS-BH except for g zz = l/z 2 + 9'(z) 2 . If 
we regard z as a "time coordinate," J = 8Cd7 /dA' x is a con- 
served "canonical momentum," since Cd7 does not contain 
A x explicitly. Then it is convenient to introduce a "Routhian" 
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C m = C m -A' x dCm/dA' x 
of A x , J, 9 and 0' full : 



lsLLTJ], which is given in terms 



9 2Z (9xx - Al/\g tt \){N*\g tt \gl x cosS 9 - J 2 ).(6) 



The Euler-Lagrange equation from © determines 9(z) under 
given E = —A x and J. The relationship between E and J 
is determined by requesting the on-shell D7-brane action to 
be real valued: the inside of the square root in ((6) has to be 
positive semidefinite iflUl . This condition yields 



NT{e 2 + 1) 



^COS 3 ! 



(7) 



where 



= [\/e 2 + 1 - e] 1/2 Zfr and e = 2£/(ttv / 2AT 2 ) 
1 1 1H : z = z* is the location where the inside of the square root 
in (0 touches zero ll 1 1. 

We need numerical analysis to obtain 6>(z») explicitly. The 
boundary conditions we employ are 9(z)/z\ z =q = m q and 
9'\ z=Zr = [B - VB 2 +C 2 ]/(Cz*). Here B = 3zfj + 
2zjjzi + 3z» and C = 3(zf - z%) tan 0(z*). The condi- 
tion for 9'\ z=z comes from the equation of motion (EOM) 
at z = z* ]5l ISD with the assumption ^ 7r/2. Note 

that #'| z=z-( is given in terms of E, T and J for given A and 
N c . After the solution 9(z) is obtained, we estimate m q from 
#(z)/z| z= o as a function of E, T and J. We choose E 1 (under 
given T and J) so that m q agrees with the designed value. 
Since the numerical analysis becomes unstable at z = 0, 
z = z*, and at z = z#, we avoid these points by introduc- 
ing small cutoffs in the numerical computations. 

Nonequilibrium phase transitions. — Let us fix N c = 40 
for numerical computation. We fix m q = 1 for simplicity. 
The J-E characteristics at various T are shown in Fig.Q] The 
system exhibits NDC in the small- J region, whereas PDC is 
seen in the large- J region. The NDC region is smoothly con- 
nected to the PDC region for T < T c , showing a crossover 
at point A. However, the curve for T > T c has an interme- 
diate region (between points D and G), where three values of 
E are possible at a given J. This multivalued nature of E at 
a given J has never been obtained in eariler works. In terms 
of the D7-brane dynamics, we have three different solutions 
to the EOM of 6(z). If we start in the small- J region, E has 
to jump to the lower value at somewhere in the intermediate 
region, and then PDC appears. Since E (hence a) changes 
discontinuously, we call it first-order transition. The bound- 
ary between the crossover regime and the first-order regime is 
found at T = 3.4365 x 10" 1 = T c . The differential resistiv- 
ity dE/dJ diverges at J = 1.86 x 10~ 2 = J c (indicated by 
B) although a changes continuously. We call this transition 
at (T, J) = (T c , J c ) a second-order transition. The minimum 
value of E for each curve is the critical electric field for insu- 
lation breaking. 

An immediate question is how to determine the transition 
point in the intermediate region. In equilibrium systems, the 
stable phase is the phase of minimum thermodynamic poten- 
tial (TP), and the transition point is where the two phases 
share a common TP. However, the generalization of the idea 
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FIG. 1: The J-E curves at T = 3.4379 x 1CT 1 > T c (circle), 
T = 3.4365 x 1CT 1 = T c (box), and T = 3.4337 x KT 1 < T c 
(triangle). Inset: the phase diagram. CP stands for the critical point. 



of TP into nonequilibrium cases has not been completely es- 
tablished. One way to evaluate the NESS generalization of TP 
is to use the Maxwell construction ll 811 . However, this method 
works only for the pairs of conjugate variables. In our case, we 
need a TP as a function of J. Note that J and E are not con- 
jugate to each other. The variable conjugate to J is A^^o, 
which explicitly depends on time, and we cannot construct a 
time-independent TP by integrating A x | z= o with respect to J: 
the method of Maxwell construction does not work for our 
purpose. Note that the Euclidean on-shell action of D7-brane 
cannot be employed as the TP in our system, since the reflec- 
tion of real time is essential in nonequilibrium systems. If we 
insist with the Euclidean formalism, we encounter unremov- 



able conical singularity in the gravity dual 019ll . In this sense, 
we do not employ the Euclidean action proposed in I20I1 . We 
need a new proposal for NESS generalization of TP within the 
Minkowski signature. 

Interestingly, we can determine the most stable phase by 
considering the Hamiltonian of the D7-brane in the gravity 
dual. The question is which D7-brane configuration is most 
stable among the three possible solutions. Since the dynamics 
of the D7-brane in the dual picture is governed by classical 
mechanics, the most stable configuration is that of the smallest 
Hamiltonian. We have seen that the D7-brane dynamics is 
governed by £07- Therefore, let us construct the Hamiltonian 
density 



%D7 = A x dC D7 /dA x - C D7 , 
which is explicitly given by I2H1 



"H D 7 = 9xx\J\gtt\gz 



'N 2 cos 6 e\g tt \g^ x - J 2 



\9tt\9x 



E 2 



(8) 



(9) 
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Note that Hd7 is regular at the horizon (where g*t = 0) 
and free from the IR divergence discussed in JII Hoj]. The 
divergence at the horizon in £ D7 is canceled by that in 
A x dCu7/dA x within the Legendre transformation. We pro- 
pose to define the thermodynamic potential per unit 3d volume 
in our system by 



H D7 (T,J;m q 



lim 



dzH 



D7 



, (10) 



where L count denotes the counterterms that renormalize the di- 
vergence at the boundary z = (which corresponds to UV 
divergence in the gauge theory). L count is given by 



■C'count — Li + L 2 — Lp + Lf, 



(11) 



where L\, Li, Lp, and Lf are explicitly given in 111 ll 122 1. 
Note that the relative sign between L p and the others has to 
be flipped, compared to the counterterms for the action owing 
to the Legendre transformation. 

It is found from our numerical analysis that the configura- 
tions between F and G have the smallest Hu7 compared to 
those between E and F, and those between D and E at a given 
J. Therefore, the transition point between NDC and PDC is D 
or F as is indicated by the arrow in Fig. [TJ Our system prefers 
the smallest dissipation under the current-controlled setup. 

The next question is how to define the critical exponents. 
We find that the NDC and PDC phases are connected via the 
crossover region on the phase diagram (shown in the inset of 
Fig-Q3 : the symmetry of the system does not change through 
the transition. This resembles the liquid-gas transitions and 
the Mott insulator-to-metal transitions in equilibrium systems, 
whose critical points are in the same universality class of the 
Ising model 12311 . In the liquid-gas transitions, the critical ex- 
ponents f3 and 5 are given by Ap oc \T — T c \° along the 
first-order transition line and \p — p c \ oc \P — Pc] 1 ^ 6 along the 
T = T c line. Here, p and P are the density and the pressure, 
and p c and P c are their critical values. Ap is the difference 
of the density between the liquid phase and the gas phase. In 
the equilibrium Mott insulator-to-metal transitions, the critical 
exponents can be detected by using the conductivity instead of 
the density |23l 12411 . Therefore, let us generalize the definition 
of the exponent /? into the nonequilibrium cases as 



A a oc \T - T c 



(12) 



where the temperature is that of the heat bath and Act is the 
difference of the conductivity between the NDC phase and the 
PDC phase along the first-order transition line. We regard Act 
as a probe of the order parameter. The pressure of the system 
is not a control parameter within the present setup, but we 
have J instead. Let us define a new critical exponent 5 by 



\a-a c \ oc \J-J C \ 1/S , 



(13) 



along the T = T C line. 

The behaviors of the conductivity are plotted in Figs. |2J a ) 
and |2jb). Critical phenomena with j3 = 0.52 ± 0.03 and 



(a) 



b 
< 



0.1 L 
10"'' 



0.0001 
(T-Tc)/Tc 



0.001 



b 
i 

b 

0.1 
0.01 

0.01 r_ 



(b) 




• • * 





0.9 



0.7 



0.3 



0.1 

(J-Jc)/Jc 



(d) 



0.001 L 
10"' 



1 

0.01 



0.001 



0.1 



(c) 




..• 


PDC ; 






Km 


Average- 


■ 




■ 






■ 

". NDCt 

■ 

■ 




■ 

■ 

■ ; 
■ 



0.0005 
(T-Tc)/Tc 



(e) 




a 





0.0001 
(T-Tc)/Tc 



0.001 0.01 



0.1 1 
(J-Jc)/Jc 



FIG. 2: Critical behaviors of various quantities, (a): Act, (b): a — a c , 
(c): conductivities and the average, (d): A{qq), and (e): (qq) c — {qq). 



5 = 3.1 ± 0.2 are numerically found. The conductivities in 
the two phases and their average, along the first-order tran- 
sition line, are shown in Fig. 12c). They resemble the coex- 
istence line and the law of rectilinear diameter in the liquid- 
gas transitions 12511 : the average conductivity shows a linear 
behavior within the numerical error. In Figs. |2}d) and 0e), 
the behaviors of the chiral condensate are shown. The chi- 
ral condensate is more sensitive to the possible numerical er- 
rors since it is read by d^9\ z= o: we need further analysis to 
estimate the precise values of the exponents. Currently, we 
observe preliminary values of the exponents /3 c hirai ~ 0.4 
and <5 c hi ra i ~ 3, where we have defined the exponents by 
A(qq) oc \T — T c \^ a6M along the first-order-transition line and 
by \{qq) - (qq)c\ <X \J - J c |V8*w a i ong the T = T c line. 
Here, A (qq) is the difference of (qq) between the two phases 
and (qq) c is the chiral condensate at the critical point. 
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